This data notebook is based on a model presented at the 2021 International Conference on Evolving Cities, University of Southampton, 22 – 24 September 2021.
If you want to cite the method/model please use:
Anderson, B. (2021). Simulating the consequences of an emissions levy at the city and neighbourhood scale. Paper presented at the International Conference on Evolving Cities, MAST Mayflower Studios, Southampton, United Kingdom. 22 - 24 Sep 2021.
If you wish to re-use material from this data notebook please cite it as:
Ben Anderson (2021) Data notebook: Simulating a local emissions levy to fund local energy effiency retrofit: Wight. University of Southampton, United Kingdom
License: CC-BY
Share, adapt, give attribution.
This data notebook estimates the value of an emissions levy using LSOA level data on emissions derived from the CREDS place-based emissions calculator. These emissions are all consumption, gas and electricity. It does this under two scenarios - a simple carbon value multiplier and a rising block tariff.
It then compares these with estimates of the cost of retrofitting EPC band dwellings D-E and F-G in each LSOA and for the whole area under study.
Key results:
This data notebook estimates a model of an emissions levy using LSOA level data on emissions derived from the CREDS place-based emissions calculator.
The model applies carbon ‘values’ to a number of emissions categories to estimate the levy revenue that would be generated for each LSOA in year 1 of such a levy. It then sums these values to given an overall levy revenue estimate for the area in the case study.
The data notebook then use estimates of the cost of retrofitting EPC band dwellings D-E and F-G together with estimates of the number of such dwellings in each of the LSOAs to calculate the likely cost of such upgrades in each LSOA and for the whole area in the case study.
Finally the data notebook compares the distributions of the two to understand whether sufficient revenue would be generated within each LSOA to enable the per-LSOA or whole case study area costs of the energy efficiency upgrades to be met. In doing so the data notebook also analyses the extent to which redistribution of revenue from high emissions areas (households) would be required.
It should be noted that this area level analysis uses mean emissions per household. It will therefore not capture within-LSOA heterogeneity in emissions and so will almost certainly underestimate the range of the household level emissions levy values that might be expected.
NB: no maps in the interests of speed
The model uses a number of Lower Layer Super Output Area (LSOA) level datasets to analyse the patterns of emissions.
All analysis is at LSOA level. Cautions on inference from area level data apply.
See https://www.creds.ac.uk/why-we-built-a-place-based-carbon-calculator/
“The highest carbon areas have an average per person footprint more than eight times larger than the lowest carbon areas.”
“We are not effectively targeting decarbonisation policies in high carbon areas. For example, the recently collapsed Green Homes Grants scheme provided a grant to cover 66% of the cost (up to £5,000) of retrofitting homes. For people claiming certain benefits, the cap was raised to 100% and £10,000. But the calculator shows that the big polluters are the large homes in very wealthy areas. In these neighbourhoods, the issue is not affordability but motivation. For high income households, energy costs are a small proportion of their expenditure and so the cost savings for retrofitting their home are inconsequential. As there are no policy “sticks” to incentivise action in the collective interest it is unsurprising that high carbon neighbourhoods have not prioritised decarbonisation."
Source: https://www.carbon.place/
Notes:
region | nLSOAs | mean_KgCo2ePerCap | sd_KgCo2ePerCap |
South East | 84 | 6,700.2 | 1,966.0 |
Now we need to convert the per capita to totals and then use the number of electricity meters as a proxy for the number of dwellings
Ideally we’d have Census 2021 data but we don’t have it yet. So instead we’ll use the number of electricity meters for 2018 which aligns with the CREDS data (might be an over-estimate if a dwelling has 2…)
First check the n electricity meters logic…
## LSOA11NM WD18NM nGasMeters nElecMeters epc_total
## 1: Isle of Wight 012H Freshwater South 1397 1361 758
## 2: Isle of Wight 001A Cowes North 1296 1729 1060
## 3: Isle of Wight 010B Brading, St Helens and Bembridge 1017 1149 564
## 4: Isle of Wight 010A Brading, St Helens and Bembridge 986 1110 606
## 5: Isle of Wight 002D Cowes West and Gurnard 986 868 632
## 6: Isle of Wight 014D Sandown South 961 1137 676
## LSOA11NM WD18NM nGasMeters nElecMeters epc_total
## 1: Isle of Wight 001A Cowes North 1296 1729 1060
## 2: Isle of Wight 012H Freshwater South 1397 1361 758
## 3: Isle of Wight 003A East Cowes 952 1203 728
## 4: Isle of Wight 010B Brading, St Helens and Bembridge 1017 1149 564
## 5: Isle of Wight 010D Nettlestone and Seaview 929 1139 604
## 6: Isle of Wight 014D Sandown South 961 1137 676
LSOA11NM | WD18NM | nGasMeters | nElecMeters | epc_total |
Isle of Wight 012H | Freshwater South | 1,397 | 1,361 | 758 |
Isle of Wight 001A | Cowes North | 1,296 | 1,729 | 1,060 |
Isle of Wight 010B | Brading, St Helens and Bembridge | 1,017 | 1,149 | 564 |
Isle of Wight 010A | Brading, St Helens and Bembridge | 986 | 1,110 | 606 |
Isle of Wight 002D | Cowes West and Gurnard | 986 | 868 | 632 |
Isle of Wight 014D | Sandown South | 961 | 1,137 | 676 |
Check that the number of electricity meters reasonably correlates with the number of EPCs from the CREDS data. We would not expect the number of gas meters to correlate due to non-gas dwellings etc.
There may also be difficulties where there are multiple meters per property - e.g. one ‘standard’ and one ‘economy 7.’ Really should switch to using address counts from postcode file.
Check that the assumption seems sensible…
Check for outliers - what might this indicate?
We want to present the analysis in ‘per dwelling’ or ‘per household’ terms so we need to convert the total kg CO2e values to per dwelling values by dividing by the number of electricity meters.
## # Summary of per dwelling values
| Name | …[] |
| Number of rows | 84 |
| Number of columns | 9 |
| Key | NULL |
| _______________________ | |
| Column type frequency: | |
| numeric | 9 |
| ________________________ | |
| Group variables | None |
Variable type: numeric
| skim_variable | n_missing | complete_rate | mean | sd | p0 | p25 | p50 | p75 | p100 | hist |
|---|---|---|---|---|---|---|---|---|---|---|
| CREDStotal_kgco2e_pdw | 0 | 1 | 13416.69 | 4470.95 | 5802.13 | 9655.92 | 12619.80 | 16100.85 | 26586.58 | ▆▇▃▂▁ |
| CREDSgas_kgco2e2018_pdw | 0 | 1 | 1924.43 | 418.34 | 817.52 | 1674.40 | 1951.28 | 2198.01 | 3018.99 | ▁▅▇▅▁ |
| CREDSelec_kgco2e2018_pdw | 0 | 1 | 1001.27 | 138.27 | 767.11 | 917.96 | 963.49 | 1047.53 | 1566.14 | ▅▇▂▁▁ |
| CREDSmeasuredHomeEnergy_kgco2e2018_pdw | 0 | 1 | 2925.70 | 444.50 | 1729.15 | 2623.09 | 2943.84 | 3236.13 | 4228.59 | ▁▇▇▅▁ |
| CREDSotherEnergy_kgco2e2011_pdw | 0 | 1 | 111.38 | 171.65 | 0.00 | 39.38 | 62.91 | 108.95 | 1141.39 | ▇▁▁▁▁ |
| CREDSallHomeEnergy_kgco2e2018_pdw | 0 | 1 | 3037.09 | 465.49 | 1800.77 | 2703.79 | 3034.90 | 3406.58 | 4337.08 | ▁▇▇▅▁ |
| CREDScar_kgco2e2018_pdw | 0 | 1 | 1668.46 | 429.29 | 786.70 | 1374.01 | 1684.29 | 1909.54 | 2765.45 | ▃▅▇▃▁ |
| CREDSvan_kgco2e2018_pdw | 0 | 1 | 270.29 | 182.75 | 63.19 | 174.86 | 225.66 | 279.58 | 1285.31 | ▇▂▁▁▁ |
| CREDSpersonalTransport_kgco2e2018_pdw | 0 | 1 | 1938.75 | 537.56 | 859.79 | 1587.01 | 1901.92 | 2144.14 | 3580.99 | ▃▇▇▂▁ |
Examine patterns of per dwelling emissions for sense.
Figure 4.1 shows the LSOA level per dwelling ‘all emissions’ in Tonnes CO2e as estimated by the CREDS tool against the Index of Multiple Deprivation (IMD) score and uses the size of the points to represent the % of dwellings with electric heating. Colour is used to represent the IMD decile where decile 1 is the 10% least deprived.
## Per dwelling T CO2e - all emissions
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
Figure 4.1: Scatter of LSOA level all consumption emissions per dwelling against IMD score
## Correlation test (pwcorr)
##
## Pearson's product-moment correlation
##
## data: selectedLsoasDT$IMDScore and selectedLsoasDT$CREDStotal_kgco2e_pdw
## t = -6.9669, df = 82, p-value = 7.435e-10
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## -0.7288856 -0.4548496
## sample estimates:
## cor
## -0.6097805
## Total emissions per dwelling (LSOA level) summary
## LSOA11CD WD18NM IMD_Decile_label All_Tco2e_per_dw
## Length:84 Length:84 4 :17 Min. : 5.802
## Class :character Class :character 5 :16 1st Qu.: 9.656
## Mode :character Mode :character 3 :15 Median :12.620
## 7 :12 Mean :13.417
## 2 : 9 3rd Qu.:16.101
## 6 : 5 Max. :26.587
## (Other):10
LSOA11CD | WD18NM | IMD_Decile_label | All_Tco2e_per_dw |
E01017334 | Parkhurst | 4 | 26.6 |
E01017323 | Arreton and Newchurch | 7 | 23.7 |
E01017296 | Arreton and Newchurch | 4 | 23.5 |
E01017301 | Cowes West and Gurnard | 9 | 23.0 |
E01017304 | Cowes Medina | 5 | 21.0 |
E01017286 | Binstead and Fishbourne | 8 | 20.9 |
LSOA11CD | WD18NM | IMD_Decile_label | All_Tco2e_per_dw |
E01017332 | Newport East | 1 (10% most deprived) | 5.8 |
E01017326 | Newport Central | 3 | 7.0 |
E01017336 | Ryde North East | 1 (10% most deprived) | 7.2 |
E01017363 | Ventnor East | 2 | 7.2 |
E01017345 | Ryde South | 2 | 8.1 |
E01017315 | Freshwater North | 3 | 8.2 |
Figure 4.2 uses the same plotting method to show emissions per dwelling due to gas use.
## Per dwelling T CO2e - gas emissions
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 817.5 1674.4 1951.3 1924.4 2198.0 3019.0
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
Figure 4.2: Scatter of LSOA level gas per dwelling emissions against IMD score
## Correlation test (pwcorr)
##
## Pearson's product-moment correlation
##
## data: selectedLsoasDT$IMDScore and selectedLsoasDT$CREDSgas_kgco2e2018_pdw
## t = -8.0335, df = 82, p-value = 5.996e-12
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## -0.7686678 -0.5237665
## sample estimates:
## cor
## -0.6636396
Figure 4.3 uses the same plotting method to show emissions per dwelling due to electricity use.
## Per dwelling T CO2e - elec emissions
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
Figure 4.3: Scatter of LSOA level elec per dwelling emissions against IMD score - who emits?
## Correlation test (pwcorr)
##
## Pearson's product-moment correlation
##
## data: selectedLsoasDT$IMDScore and selectedLsoasDT$CREDSelec_kgco2e2018_pdw
## t = -3.0053, df = 82, p-value = 0.003519
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## -0.4958962 -0.1078802
## sample estimates:
## cor
## -0.31499
Figure 4.4 uses the same plotting method to show emissions per dwelling due to other energy use. This should be higher for off-gas areas which tend to be rural areas so we also present analysis by the LSOA’s urban/rural classification.
## Per dwelling T CO2e - elec emissions
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
Figure 4.4: Scatter of LSOA level other energy per dwelling emissions against IMD score - who emits?
## Correlation test (pwcorr)
##
## Pearson's product-moment correlation
##
## data: selectedLsoasDT$IMDScore and selectedLsoasDT$CREDSelec_kgco2e2018_pdw
## t = -3.0053, df = 82, p-value = 0.003519
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## -0.4958962 -0.1078802
## sample estimates:
## cor
## -0.31499
RUC11 | mean_gas_kgco2e | mean_elec_kgco2e | mean_other_energy_kgco2e |
Rural town and fringe | 2,028.1 | 1,025.6 | 109.4 |
Rural village and dispersed | 1,590.2 | 1,319.2 | 597.7 |
Urban city and town | 1,921.2 | 968.0 | 72.1 |
Check whether all measured energy emissions combined (gas & electricity) correlate with all emissions (in this data).
##
## Pearson's product-moment correlation
##
## data: selectedLsoasDT$CREDStotal_kgco2e_pdw and selectedLsoasDT$CREDSmeasuredHomeEnergy_kgco2e2018_pdw
## t = 6.0727, df = 82, p-value = 3.752e-08
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## 0.3890300 0.6890812
## sample estimates:
## cor
## 0.5569704
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
Do we see strong correlations? If so in theory we could (currently) use measured energy emissions as a proxy for total emissions.
Repeat for all home energy - includes estimates of emissions from oil etc
##
## Pearson's product-moment correlation
##
## data: selectedLsoasDT$CREDStotal_kgco2e_pdw and selectedLsoasDT$CREDSallHomeEnergy_kgco2e2018_pdw
## t = 8.4743, df = 82, p-value = 7.996e-13
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## 0.5493802 0.7829857
## sample estimates:
## cor
## 0.683294
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
How does the correlation look now?
We don’t expect to use this data as it is already taxed in a way that relates to emissions (?)
Figure 4.5 uses the same plotting method to show emissions per dwelling due to van use. Again, we present analysis by the LSOA’s urban/rural classification.
## Per dwelling T CO2e - car emissions
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
Figure 4.5: Scatter of LSOA level car use per dwelling emissions against IMD score
## Correlation with IMD score (pwcorr)
##
## Pearson's product-moment correlation
##
## data: selectedLsoasDT$IMDScore and selectedLsoasDT$CREDScar_kgco2e2018_pdw
## t = -6.8002, df = 82, p-value = 1.561e-09
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## -0.7219420 -0.4431483
## sample estimates:
## cor
## -0.6004914
RUC11 | mean_car_kgco2e | mean_van_kgco2e |
Rural town and fringe | 1,799.4 | 232.2 |
Rural village and dispersed | 2,447.3 | 524.0 |
Urban city and town | 1,566.0 | 260.7 |
Figure 4.6 uses the same plotting method to show emissions per dwelling due to van use.
## Per dwelling T CO2e - van emissions
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
Figure 4.6: Scatter of LSOA level van use per dwelling emissions against IMD score
## Correlation with IMD score (pwcorr)
##
## Pearson's product-moment correlation
##
## data: selectedLsoasDT$IMDScore and selectedLsoasDT$CREDSvan_kgco2e2018_pdw
## t = -1.5383, df = 82, p-value = 0.1278
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## -0.36863496 0.04866661
## sample estimates:
## cor
## -0.167476
Case studies:
BEIS/ETC Carbon ‘price’
EU carbon ‘price’
Scenario 1: apply the central value Scenario 2: apply the low/central/high as a rising block tariff for each emissions source. Set threhsolds to 33% and 66% (in absence of any other guidance!)
Table 4.6 below shows the overall £ GBP total for the case study area in £M under Scenario 1.
nLSOAs | beis_GBPtotal_c | beis_total_c_gas | beis_GBPtotal_c_elec |
84 | 216.5 | 32.3 | 16.7 |
region | nLSOAs | beis_GBPtotal_c | beis_total_c_gas | beis_GBPtotal_c_elec |
South East | 84 | 216.5 | 32.3 | 16.7 |
Figure 4.7: Proportion of total emissions due to gas & electricity use by region covered
The table below shows the mean per dwelling value rounded to the nearest £10.
All_emissions | Gas | Electricity | Gas + Electricity |
3,287.1 | 471.5 | 245.3 | 716.8 |
Figure ?? shows the total £k per LSOA and £ per dwelling revenue using BEIS central carbon price plotted against IMD score. The tables show the LSOAs with the highest and lowest values.
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
Figure 4.8: £k per LSOA revenue using BEIS central carbon price
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
Figure 4.9: £k per LSOA revenue using BEIS central carbon price
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 1422 2366 3092 3287 3945 6514
LSOA11CD | LSOA01NM | WD18NM | T_CO2e_pdw | GBP_allEmissions_levy |
E01017334 | Isle of Wight 008E | Parkhurst | 26.6 | 6,513.7 |
E01017323 | Isle of Wight 015C | Arreton and Newchurch | 23.7 | 5,800.5 |
E01017296 | Isle of Wight 011C | Arreton and Newchurch | 23.5 | 5,763.6 |
E01017301 | Isle of Wight 002A | Cowes West and Gurnard | 23.0 | 5,626.4 |
E01017304 | Isle of Wight 002C | Cowes Medina | 21.0 | 5,154.8 |
E01017286 | Isle of Wight 005A | Binstead and Fishbourne | 20.9 | 5,128.3 |
LSOA11CD | LSOA01NM | WD18NM | T_CO2e_pdw | GBP_allEmissions_levy |
E01017332 | Isle of Wight 008D | Newport East | 5.8 | 1,421.5 |
E01017326 | Isle of Wight 009C | Newport Central | 7.0 | 1,719.2 |
E01017336 | Isle of Wight 004B | Ryde North East | 7.2 | 1,767.9 |
E01017363 | Isle of Wight 018A | Ventnor East | 7.2 | 1,769.1 |
E01017345 | Isle of Wight 007E | Ryde South | 8.1 | 1,976.4 |
E01017315 | Isle of Wight 012D | Freshwater North | 8.2 | 2,012.7 |
Figure ?? repeats the analysis but just for gas.
Anything unusual?
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
Figure 4.10: £k per LSOA incurred via gas using BEIS central carbon price
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
Figure 4.11: £k per LSOA incurred via gas using BEIS central carbon price
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 200.3 410.2 478.1 471.5 538.5 739.7
LSOA11CD | LSOA01NM | WD18NM | gas_T_CO2e_pdw | GBP_gas_levy_perdw |
E01017286 | Isle of Wight 005A | Binstead and Fishbourne | 3.0 | 739.7 |
E01017316 | Isle of Wight 002D | Cowes West and Gurnard | 2.7 | 670.0 |
E01017352 | Isle of Wight 010E | Nettlestone and Seaview | 2.7 | 657.0 |
E01017284 | Isle of Wight 010A | Brading, St Helens and Bembridge | 2.6 | 626.6 |
E01017323 | Isle of Wight 015C | Arreton and Newchurch | 2.5 | 623.9 |
E01033239 | Isle of Wight 012H | Freshwater South | 2.5 | 615.7 |
LSOA11CD | LSOA01NM | WD18NM | gas_T_CO2e_pdw | GBP_gas_levy_perdw |
E01033237 | Isle of Wight 003E | Whippingham and Osborne | 1.3 | 330.7 |
E01017325 | Isle of Wight 009B | Newport Central | 1.3 | 315.1 |
E01017345 | Isle of Wight 007E | Ryde South | 1.2 | 301.5 |
E01017332 | Isle of Wight 008D | Newport East | 0.9 | 221.0 |
E01017296 | Isle of Wight 011C | Arreton and Newchurch | 0.8 | 201.2 |
E01017353 | Isle of Wight 013D | West Wight | 0.8 | 200.3 |
Figure ?? repeats the analysis for electricity.
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
Figure 4.12: £k per LSOA incurred via electricity using BEIS central carbon price
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
Figure 4.13: £k per LSOA incurred via electricity using BEIS central carbon price
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 187.9 224.9 236.1 245.3 256.6 383.7
LSOA11CD | LSOA01NM | WD18NM | elec_T_CO2e_pdw | GBP_elec_levy_perdw |
E01017296 | Isle of Wight 011C | Arreton and Newchurch | 1.6 | 383.7 |
E01017353 | Isle of Wight 013D | West Wight | 1.5 | 364.2 |
E01017369 | Isle of Wight 017D | Godshill and Wroxall | 1.4 | 331.9 |
E01017365 | Isle of Wight 018C | Ventnor West | 1.3 | 313.8 |
E01017310 | Isle of Wight 008A | Newport North | 1.2 | 305.2 |
E01017323 | Isle of Wight 015C | Arreton and Newchurch | 1.2 | 302.7 |
LSOA11CD | LSOA01NM | WD18NM | elec_T_CO2e_pdw | GBP_elec_levy_perdw |
E01017363 | Isle of Wight 018A | Ventnor East | 0.8 | 207.5 |
E01033237 | Isle of Wight 003E | Whippingham and Osborne | 0.8 | 206.1 |
E01017303 | Isle of Wight 001B | Cowes North | 0.8 | 204.3 |
E01017305 | Isle of Wight 001C | Cowes Medina | 0.8 | 203.6 |
E01017332 | Isle of Wight 008D | Newport East | 0.8 | 202.7 |
E01017320 | Isle of Wight 015B | Lake South | 0.8 | 187.9 |
Figure ?? shows the same analysis for measured energy (elec + gas)
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
Figure 4.14: £k per LSOA incurred via electricity and gas using BEIS central carbon price
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
Figure 4.15: £k per LSOA incurred via electricity and gas using BEIS central carbon price
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 423.6 642.7 721.2 716.8 792.9 1036.0
Applied to per dwelling values (not LSOA total) - may be methodologically dubious?
Cut at 25%, 50% - so any emissions over 50% get high carbon cost
## Cuts for total per dw
## 0% 25% 50% 75% 100%
## 5802.128 9655.919 12619.805 16100.849 26586.576
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.
Table: (#tab:estimateAnnualLevyScenario2Total)Data summary
| Name | …[] |
| Number of rows | 84 |
| Number of columns | 3 |
| Key | NULL |
| _______________________ | |
| Column type frequency: | |
| numeric | 3 |
| ________________________ | |
| Group variables | None |
Variable type: numeric
| skim_variable | n_missing | complete_rate | mean | sd | p0 | p25 | p50 | p75 | p100 | hist |
|---|---|---|---|---|---|---|---|---|---|---|
| V1 | 0 | 1 | 13.42 | 4.47 | 5.80 | 9.66 | 12.62 | 16.10 | 26.59 | ▆▇▃▂▁ |
| beis_GBPtotal_sc2_perdw | 0 | 1 | 2392.12 | 1447.48 | 707.86 | 1179.62 | 1906.62 | 3181.72 | 7029.98 | ▇▃▂▂▁ |
| beis_GBPtotal_sc2 | 0 | 1 | 1813028.61 | 868870.61 | 665388.00 | 1098649.80 | 1617217.51 | 2481316.82 | 3914149.60 | ▇▅▃▃▂ |
nLSOAs | sum_total_sc1 | sum_total_sc2 |
84 | 216.5 | 152.3 |
Figure 4.16 compares the % £ levy under each scenario for all consumption contributed by LSOAs in each IMD decile.
Figure 4.16: Comparing £ levy under each scenario by IMD decile - all consumption emissions
Figure 4.17 compares the £ levy under each scenario for all consumption.
Figure 4.17: Comparing £ levy under each scenario - all consumption emissions
## [1] 20.00895
## [1] 10.26845
nLSOAs | sumAllConsEmissions_GBP | sumGasEmissions_GBP | sumElecEmissions_GBP |
84 | 152.3 | 20.0 | 10.3 |
region | nLSOAs | sumAllConsEmissions_GBP | sumGasEmissions_GBP | sumElecEmissions_GBP | sumPop |
South East | 84 | 152.3 | 20.0 | 10.3 | 134,820 |
Figure 4.18: Contribution to sum levy £ GBP by source
Source: English Housing Survey 2018 Energy Report
Model excludes EPC A, B & C (assumes no need to upgrade)
Adding these back in would increase the cost… obvs
In order to estimate the LSOA level retrofit costs, we need to impute the EPC counts in each LSOA. We do this using the number of electricity meters as the presumed number of dwellings and the observed % of EPCs in each band for all dwellings with EPCs which is provided by the CREDS data. This assumes that if we had EPCs for all dwellings then the % in each band in each LSOA would stay the same. This is quite a bold assumption…
Note that the EPC database is continuously updated so more recent upgrades will not be captured in the data used for this analysis. This means the total retrofit costs are likely to be an over-estimate. The extent of this over-estimate would require the use of an updated (current) EPC data extract and is left for future work.
## N EPCs
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 231.0 372.0 421.0 464.7 549.0 1060.0
## N elec meters
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 513.0 692.5 773.0 817.4 896.2 1729.0
Correlation between high % EPC F/G or A/B and deprivation?
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
Now we need to convert the % to dwellings using the number of electricity meters (see above).
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
Table 4.13 reports total retofit costs.
## To retrofit D-E (£m)
## [1] 504.9341
## Number of dwellings: 37965
## To retrofit F-G (£m)
## [1] 126.1724
## Number of dwellings: 4708
## To retrofit D-G (£m)
## [1] 631.1064
## To retrofit D-G (mean per dwelling)
## [1] 14689.25
meanPerLSOA_GBPm | total_GBPm |
7.5 | 631.1 |
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
Figure 4.19 shows the LSOA level retofit costs per dwelling by IMD decile.
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
Figure 4.19: LSOA level retofit costs per dwelling by IMD score
Totals
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
Repeat per dwelling
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
Totals
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
Repeat per dwelling
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
Figure 4.20 shows years to pay under Scenario 1 (all emissions)
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 2.197 3.718 4.710 4.959 6.251 9.923
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
Figure 4.20: Years to pay under Scenario 1 (all em issions)
## Median years: 4.71
Figure 4.21 shows years to pay under Scenario 1 (energy emissions)
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 14.05 18.55 20.54 20.99 23.05 33.30
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
Figure 4.21: Years to pay under Scenario 1 (energy emissions)
## Median years: 20.54
Figure 4.22 shows the year 1 outcome if levy is shared equally (all emissions levy).
Figure 4.22: Year 1 outcome if levy is shared equally (all emissions levy)
LSOA11CD | LSOA11NM | WD18NM | retrofitSum | yearsToPay | epc_D_pc | epc_E_pc | epc_F_pc | epc_G_pc |
E01017300 | Isle of Wight 001A | Cowes North | 15,404,574.4 | 22.0 | 0.4 | 0.2 | 0.1 | 0.0 |
E01033239 | Isle of Wight 012H | Freshwater South | 14,972,795.5 | 17.3 | 0.4 | 0.2 | 0.1 | 0.0 |
E01017351 | Isle of Wight 010D | Nettlestone and Seaview | 14,534,696.0 | 20.6 | 0.3 | 0.3 | 0.1 | 0.0 |
E01017337 | Isle of Wight 004C | Ryde North West | 12,228,840.0 | 24.0 | 0.4 | 0.3 | 0.1 | 0.0 |
E01017365 | Isle of Wight 018C | Ventnor West | 12,211,760.7 | 24.5 | 0.3 | 0.2 | 0.1 | 0.1 |
E01017363 | Isle of Wight 018A | Ventnor East | 11,758,862.1 | 28.6 | 0.3 | 0.2 | 0.1 | 0.0 |
E01017367 | Isle of Wight 005C | Wootton Bridge | 11,531,191.8 | 18.5 | 0.5 | 0.2 | 0.1 | 0.0 |
E01017364 | Isle of Wight 018B | Ventnor East | 11,465,197.9 | 20.5 | 0.4 | 0.3 | 0.1 | 0.1 |
E01017336 | Isle of Wight 004B | Ryde North East | 11,417,419.0 | 28.2 | 0.3 | 0.2 | 0.1 | 0.0 |
E01017285 | Isle of Wight 010B | Brading, St Helens and Bembridge | 10,935,261.2 | 16.9 | 0.4 | 0.2 | 0.1 | 0.0 |
Figure 4.23 shows the year 1 outcome if levy is shared equally (energy emissions levy).
Figure 4.23: Year 1 outcome if levy is shared equally (energy emissions levy)
What happens in Year 2 totally depends on the rate of upgrades… given the supply chain & capacity issues it’s likely that the levy would build up a substantial ‘headroom’ that could then be spent over time…
Figure 4.24 shows years to pay under Scenario 2 (all emissions)
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 2.036 4.635 7.741 8.357 12.472 19.928
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
Figure 4.24: Years to pay under Scenario 2 (all em issions)
## Median years: 7.74
Figure 4.25 shows years to pay under Scenario 2 (energy emissions)
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 14.05 18.55 20.54 20.99 23.05 33.30
## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
Figure 4.25: Years to pay under Scenario 2 (energy emissions)
Figure 4.26 shows the year 1 outcome if levy is shared equally (all emissions levy).
Figure 4.26: Year 1 outcome if levy is shared equally (all emissions levy)
Figure 4.27 shows the year 1 outcome if levy is shared equally (energy emissions levy).
Figure 4.27: Year 1 outcome if levy is shared equally (energy emissions levy)
What happens in Year 2 totally depends on the rate of upgrades…
Figure 4.28 compares pay-back times for the two scenarios - who does the rising block tariff help?
Figure 4.28: Comparing pay-back times across scenarios
I don’t know if this will work…
## Doesn't